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SUMMARY 


The development of shocks in the sound field propagating 
through a nearly choked duct flow is analyzed by extending the 
quasi-one dimensional nonlinear theory previously derived by the 
authors. The theory is applied to the case in which sound is in- 
troduced into the flow by an acoustic source located in the vicin- 
ity of a near-sonic throat. Analytical solutions for the field 
are obtained which illustrate the essential features of the non- 
linear interaction between sound and flow. 

Numerical results are presented covering ranges of variation 
of source strength, throat Mach number, and frequency. It is 
found that the development of shocks leads to appreciable attenua- 
tion of acoustic power transmitted upstream through the near-sonic 
flow. It is possible, for example, that the power loss in the 
fundamental harmonic can be as much as 90% of that introduced at 
the source. 


1. INTRODUCTION 


The possibility of reducing noise transmitted upstream in 
engine inlets by using a nearly choked mean flow has been recognized 
for some time. In a sequence of papers (refs. 1-6) , the present 
authors have developed a nonlinear quasi-one dimensional theory of 
sound propagation through a high-subsonic throat which is aimed at 
understanding the physical mechanisms responsible for the noise 
reduction. The theory is derived using the method of matched 
asymptotic expansions. With the exception of the work reported in 
reference 6, the theory has thus far been applied to cases in which 
shock waves do not occur. It is clear, however, that it offers a 
possible explanation for a reduction in upstream travelling sound. 

A near-sonic steady flow near the throat causes a transonic trap- 
ping of waves propagating upstream, and the resulting intensifica- 
tion of the field gives rise to nonlinear interactions between the 
sound and the flow. As a result, for example, sound energy intro- 
duced by a harmonic source is partitioned into an infinite number 
of superharmonic components and a steady streaming component after 
passing through the near-sonic flow region. Even in the absence 
of any actual dissipation of energy, the energy transmitted out of 
the throat region in the fundamental harmonic can be reduced sig- 
nificantly compared to that imparted at the source. 

Previous work by other authors which relates to the present 
investigation has been discussed in references 1-6. No attempt 
will be made here to survey all of these sources again. However, 
it should be noted that a fully numerical investigation of the non- 
linear quasi-one dimensional problem has been carried out by Nayfeh 
and his co-workers (see refs. 7,8). Their numerical results con- 
firm the existence of all of the nonlinear effects mentioned above 
in the shock-free case. 

By far the most significant physical consequence of the non- 
linearity, however, is the fact that sufficiently high source 
strengths or frequencies cause the nonlinear acoustic field to 
develop shock waves in the near-sonic flow region which ultimately 
appear in the sound field away from the throat region as simple 
acoustic discontinuities in each period of the sound wave. In 
earlier work (refs. 1-5) this effect was recognized, and the condi- 
tions under which shocks occur were derived. It is the purpose of 
the current report to carry out the necessary generalization of the 
theory to account for these shock waves. As will be seen, this 
involves an extensive nonlinear analysis which shows that compli- 
cated shock motions can occur in the throat region. The shocks are 
associated with two related physical effects on the sound field. 
First, they actually dissipate sound energy through the thermo- 
viscous effects modelled by the shock jump conditions. Second, as 
will be seen in detail later, sound introduced during relatively 
high amplitude portions of the period can interact with the flow to 
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produce local supersonic conditions at the source during other 
portions of the period. The flow can then be said to be acous- 
tically choked: the sound signal is blocked during this portion 

of the period from propagating upstream. It will be seen in the 
following that these nonlinear effects indeed produce a signifi- 
cant reduction of sound transmitted out of the near-sonic flow 
region. At present, to the authors' knowledge, there has been no 
previous analytical or numerical treatment of the shocked acoustic 
field reported in the literature. 

The analysis in the present work is concerned primarily with 
the situation in which sound is introduced into the flow in the 
vicinity of the throat of a converging-diverging duct. It is a 
generalization of the work which appears in reference 6. This 
arrangement illustrates essentially all of the nonlinear effects 
of a near-sonic flow on sound, and it offers the advantage over a 
more general case that the solution for the field can be obtained 
analytically. Numerical results are given as an illustration of 
the theory for the special case in which the sound source is loca- 
ted exactly at the throat of the duct. 

In Chapters 3 and 4 of the report the analytical solution for 
the linearized acoustic field away from the throat is reviewed, 
and the basic nonlinear theory valid near the throat is set out in 
some detail. The solution to the nonlinear inner problem is de- 
veloped in Chapter 5, including a detailed analysis of the shock 
waves which arise in the throat region. Then in Chapter 6 the 
shocked inner solution is asymptotically matched to the outer sol- 
ution to yield a complete description of the field in terms of a 
set of easily evaluated integrals and the solution of two pairs of 
transcendental equations. 

One major result of the present analysis is the development 
of a measure of the sound energy dissipated by shocks in the 
throat region of the duct. The usual concept of acoustic energy 
flux is not sufficient to study sound energy in this region owing 
to the fact that the field there is not described by linearized 
theory. In Chapter 7 an appropriate generalization of acoustic 
energy flux is introduced, and this leads to a relatively simple 
formula describing the sound energy dissipated by thermoviscous 
effects in the nonlinear theory. This dissipation process is ex- 
tremely significant: typical results presented later show that 

as much as 40 to 50% of the total sound power introduced at the 
source can be lost during propagation of the sound upstream when 
the throat Mach number is only as high as 0.9. In combination 
with this dissipation, the spreading of energy among harmonics can 
lead to a reduction of power in the fundamental in the vicinity of 
90%. It is emphasized that this almost total shutting down of 
sound in the first harmonic occurs even though the basic steady 
flow is relatively far from being choked. 
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2 . 


SYMBOLS AND NOMENCLATURE 


A(x) 

A± ,b£ 
In' In 


n 


c 0 (x) 

C 11 ' C 12 ' d 02 
C n (x) ,D n (x) 

D(N) 

E , E* , E , E 
g (x; e) 

G ( x ; e ) 

G o (x) 

G^ n , i-1 , • • . , 4 
H (x; e) 

V x) 

I 

I, ,I 0 ,J 
In 2n n 

K 


L 

M ( x ; e ) 

M q (x) 
m 1 (X) 

N , N 

p ( x , t ; e ) 
p(x, t ; e) 

Pl (x;t; e) 

P 0 (x,t) 

r x (X , t) ,r 2 (X f t) 
R(x) 


throat curvature parameter 

2 

duct cross-sectional area, m 
see equation (3.11) 
magnitude of A^ n 
stagnation sound speed, ms ^ 
sound speed in steady flow 
see equation (3.12) 

Fourier coefficients (eq. 6.4) 
integration constant (eq. 5.24) 
dimensionless energy densities 
see equation (6.3) 
see equation (3.5) 

G (x; e=0) 

see equation (3.12) 
shock location 

leading order inner representation of H(x;£) 
time-averaged dimensionless power flow 
see equations (6.9-6.10) 

parameter in Crocco-Tsien Mach number 
distribution 

characteristic length, m 
steady flow Mach number 
M (x; e=0) 

inner expansion Mach number term 
source amplitudes 

-2 

total fluid pressure, Nm 
dimensionless perturbation pressure 
first-order outer perturbation pressure 
P x (x, t ? e=0) 

inner density expansion terms 

-3 

steady flow density, kgm 
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dimensionless time 


t 

T 

U,U,U 
u(x,t;e) 
u(x; t; e) 
u 1 ( x , t ; £ ) 
u Q (x, t) 

V 

W,W*,W,W 

"3 

Wi 

n 

x 


x 0 ,x e 


X 

X 

Y 

6 


£ 

e 

n 

y, (X, t) ,y 9 (X,t) 

n(x,t) 

p ( x , t ; e ) 

P (X,t?£) 
pj( x , t ; e ) 

P Q (x, t) 

o 1 (x,t) ,a 2 (x,t) 

T 

<MT) 

Y n 

0) 


modified dimensionless time, a^t 
source velocity amplitude ms 1 
total fluid velocity, ms~^ 
dimensionless perturbation velocity 
first order outer perturbation velocity 

u, (x, t ; e=0) 

x _1 
shock velocity, ms 

dimensionless energy fluxes 

first order outer representation of W 

inner expansion energy flux terms 

dimensionless axial coordinate 

source and exit locations 

inner axial coordinate (eq. 4.1) 

inner location of source 

ratio of specific heats 

outer acoustic disturbance level 

throat Mach number deviation from unity 

phase angle in inner solution 

inner velocity expansion terms 

inner dependent variable (eq. 5.3) 

_ 3 

total fluid density kgm 
dimensionless perturbation density 
first-order outer perturbation density 
P 1 (x, t;£=0) 

inner pressure expansion terms 

characteristic parameter 

source function (eq. 5.10) 

phase angle in outer expansion (eq. 6.1) 

dimensionless frequency 

-J- 

modified frequency, a 2 u) 
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Special Notation 


< > 

[ ] 

subscript s 
subscript n 
superscript i 
superscript ° 
superscript + 
subscript 1,2 


denotes magnitude 

denotes time average over one period 

denotes jump across shock 

denotes stagnation value 

denotes a harmonic 

denotes inner variable 

denotes outer variable 

denotes x>0, x<0, respectively 

used on t,<J>(t), IT (t ) to denote downstream 
and upstream sides of shock, respectively 
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3. FORMULATION AND OUTER EXPANSIONS 


Perturbation Equations 

Provided that the acoustic wavelength is sufficiently large 
and the area variation is sufficiently slow, the propagation of 
sound in a duct of cross-sectional area A(x) carrying an inviscid 
ideal gas flow can be described by the equations of quasi-one 
dimensional gas dynamics (ref. 9) : 


C s^t 

+ 

^x 

+ 

PU x + pu A'/A = o. 

c s“t 

+ 

-X 

+ 

P x /p = 0, 

C sP t 

+ 

“P X 

- 

c 2 (c s P t + up x ) = 0 . 


In equation (3.1) P, u, and p are the total fluid density, axial 
velocity, and pressure, respectively; c is the sound speed in the 
gas, and c s is its stagnation value. The independent variables x 
and t are dimensionless, being measured in units of L and L/c s , 
respectively, where L is a characteristic length associated with 
the axial variation in A(x). The geometry of the problem is as 
indicated in figure 1 ; the origin of x is chosen to correspond to 
a throat: A' (0) = 0. 


The velocity, density, and sound speed of the basic steady 
flow are denoted by U(x), R(x) , and cg(x), respectively. Then the 
pressure in the steady flow is cq^R/y, where y is the specific 
heat ratio of the gas, and U(x) and R(x) satisfy 


UR' + RU' + RUA'/A = 0 
UU' + c 0 2 R'/R = 0. 


(3.2) 


It is convenient to study unsteady perturbations about the 
basic steady flow in the duct in terms of dimensionless perturba- 
tion quantities p(x,t), u(x,t), and p(x,t) defined by 

p(x,t) = R (x) [ 1 + p (x , t) ] , u ( x , t ) = U (x) [ 1 + u(x,t)], 

c 0 2 r (3.3) 

p (x , t) = — — [1 + p(x,t) ] . 

The usual acoustic perturbations are then given by R(x)p(x,t), 
U(x)u(x,t), and cq^(x)R(x) p(x,t)/y. Substitution of equations 
(3.3) into equations (3.1) and employing the relations (3.2) yields 
the system of equations on p, u, and p in the form 
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G^P t + M{(l+u)p x + (l+p)u x > = 0, 

G '% + M(l+u)u x + ymIT+pT” Px + Ir^d+u) 2 - - »- 

G% t + M(l+u)p x - Y .[j:pP ) {G h P t + M(1+U)P X > = 0, 


(3.4) 


in which M(x) is the Mach number of the basic steady flow, 
U(x)/Cq(x), and 

G(x) = c s 2 /c 0 2 (x) = 1 + (y-l)M 2 /2. (3.5) 

It will also be required to make use of the shock relations 
appropriate to discontinuous solutions p, u, and p. These are 
(ref. 10) , 

-V[p] + [pu] = 0, 

-V[pu] + [pu 2 + p] = 0, (3.6) 

- V[ ^r + ylTPl + + ft^] = °' 

in which V is the velocity of the shock , and the square brackets 
denote the jump in the enclosed quantity across the shock. Sub- 
stituting the definitions (3.3) into conditions (3.6) and assum- 
ing that the basic steady flow is continuous leads to the jump 
conditions on the perturbation quantities in the form 

-V[p] + U[p + u + pu] = 0, 

+ u 2 ] + ... =0, 

(3.7) 


0 , 


and energy jump conditions 
are third-order products of the perturbation quantities which are 
not needed in the following. 


+ u + pu] + u[p + 2u + + 2 pu 

yM 

+ U + 2 + pu + ^j-] 

y(y-l)M^ 

+ U[(| + i — ~)u + §• + 2 — 

* (y-l)M z z (y-l)M z 

+ pu + — ~-] + . . . = 

(y-l)M 

where the terms omitted in the momentum 


-V[p 


-V[f 


It is assumed that the basic flow quantities R, U, and c Q 
are known from solving the system (3.2) subject to suitable 
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boundary conditions applied far upstream and downstream of the 
duct throat at x=0 . Of interest in the present work is the solu- 
tion of the quasi-linear system (3.4), which describes the per- 
turbation quantities exactly, regardless of their magnitude. In 
general these quantities will contain discontinuities whose mag- 
nitude and space-time path are determined from the shock condi- 
tions (3.7) . 

Because sound propagation in near-sonic flows is of interest 
here, it will be assumed that the boundary conditions on the basic 
steady flow are such that the throat Mach number is near unity. 

A parameter e=l-M(0) is thus introduced into the analysis, and the 
Mach number is henceforth regarded as a function of the parameter 
e as well as of x. In the analysis that follows e will be assumed 
to be positive: the basic flow is subsonic throughout the duct. 

The coefficients in the system (3.4) and in the conditions (3.7) 
then all become functions of e, which implies that the perturba- 
tion quantities are functions of e. In the work which follows, 
asymptotic solutions for p, u, and p will be derived in the limit 
as £-*0. It has been shown in detail in previous publications by 
the authors (refs. 2-5) that for e<<l the acoustic perturbation 
field can be divided into two regions: an outer region away from 

the throat where linear acoustic theory is a valid first approxi- 
mation to the system (3.4), and an inner region around the throat 
where nonlinear effects predominate. The essential results of the 
previous analyses of the behavior of the acoustic quantities in 
the outer region will be summarized in the following section. 


Outer Expansions: Linearized Acoustic Theory 


If it is assumed that the perturbation quantities p, u and p 
are small, the system (3.4) can be linearized. Let 5 be a small 
parameter measuring the order of magnitude of the perturbations. 
Expressing the perturbation quantities in powers of 5 as 


p (x,t; e , 6) = 6p 1 (x,t;£) + 0(6 2 ), 
u(x,t;e,6) = 6u^(x,t;£) + 0(6 2 ), 
p(x,t?£,<5) = 6p 1 (x,t;£) + 0 (S 2 ) , 


(3.8) 


substituting equations (3.8) into equations (3.4), and neglecting 
all but the terms proportional to 6 yields the linearized acous- 
tic equations 


10 



(3.9) 


gS, + M(p. + u. ) = 0, 
x t x x x 

G V + Mu i + w p i + Sr (2u i + p i ' p i> - °' 

t X X 

G^p + Mp. - T(G Js P, + Mp n ) =0. 

t X X X t X X 

In general, the system (3.9) would have to be solved numer- 
ically, subject to appropriate boundary and initial conditions, 
because of its variable coefficients. For the present purposes, 
however, since £<<1, a further expansion of the perturbation 
quantities is carried out such that 

P 1 (x,t;e) = p 0 (x , t ) + 0(e), 

u^ ( x , t ; e ) = Uq (x, t) + 0(e), 

P 1 (x,t;e) = p Q (x,t) + 0(e), 

where the quantities pg, Uq, and Pq satisfy the linear system 
(3.9) with the coefficient terms M(x;e) and G(x;e) replaced by 

Mq (x) = M(x;0) and Gq(x) = G(x;0). 

As has been discussed previously (ref. l), the function Mq(x) is 
any Mach number distribution such that Mg(0)=l and Mg(x)<l 
throughout the remainder of the duct. In general, M Q (x) has the 
expansion for x<< 1 (ref. 1) 

Mg (x) = 1 - a}' 2 | x | + 0 (x 2 ) (3.10) 

in which a = A" (0) /2A (0) >0 . 

The solution for Pg, ug and pg for a general area distribu- 
tion has been discussed extensively by the authors (refs. 1,5). For 
the present work, however, a specific duct shape will be chosen 
for which the authors obtained an analytical solution by extending 
a result of Crocco and Tsien (ref. 1) . The Mach number distribu- 
tion for the Crocco-Tsien duct is given for all x by 

Mg (x) = (1-K|x| ) (1-K|x|) 2 f1 | x [ <~. 

Near x=0 this assumes the form of equation (3.10) with 

= 2a/ (y+1) . If a dimensionless frequency to, measured in units 
of c s /L is introduced, then complex periodic solutions to the sys- 
tem (3.9) to leading order in e are 
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oo 


U , 


= E 
n=0 


= E 
n=0 


{A ln G ln (x > 


{fl ln G 3n (x) 


+ B ln G 2n (x > ,e 


inwt 


(3.11) 


B ln G 4n (x),e 


inwt 


P 0 = YP 0' 

where A^ n and B^ n are arbitrary complex constants, and the G^ n , 

i=l,..., 4, are combinations of hypergeometric functions of 

complex argument whose explicit forms are given in the appendix. 

The local behavior of the G. as x+0 has been shown to be (refs. 1,5), 

m 


(■pvT + c n ) exp sgn x ln|x|) 

a 2 


j-, (-1 r + C-, -| J 

In | x J 11 


Sn ^ 1 


’3n ~ 1 fxT C 12 1 


(“T5T + c i^ ex P (^jr- sgn x ln|x|) 

a 2 


(3.12) 


G 4n ~ d 02 


where the constants c.^, c±2' an< ^ ^02 a PP ear also in the appendix. 

The local expansion (3.12) provides the key to the analysis 
of the acoustic behavior of the converging-diverging duct. The 
acoustic quantities are given by 

p = 6p 0 ( x , t ) + 0 (eS) 

u = 6 uq ( x, t) + O (eS) (3.13) 

p = 5y Pq ( x , t ) + 0(e6) 

in the limit as e->0 . For any 6, no matter how small, the linear- 
ized theory fails to approximate the perturbation quantities near 
the throat as the Mach number there approaches unity. Equations 
(3.12) show that the leading terms p Q , Ug , and p 0 of the expan- 
sions (3.13) become arbitrarily large in this circumstance, there- 
by violating the assumption made in deriving the system (3.9) that 
they remain small. Numerical solutions of the system (3.9) which 
illustrate the growth of these large amplitudes were presented in 
reference 1. 
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Thus the interpretation of the role of linearized acoustic 
theory in the current problem is clear. The expansion (3.8) is 
to be regarded as an outer asymptotic expansion, valid as 6+0 with 
x, t, and e fixed. The result (3.12) shows that it is not uni- 
formly valid as x-*0 and £->-0 . Hence, in the vicinity of the duct 
throat, nonlinear effects become important as the flow approaches 
sonic conditions. The correct nonlinear theory which must be 
used to describe sound in the near-sonic throat flow will be pre- 
sented in the next chapter. 
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4 . NONLINEAR THEORY 


Inner Expansions 


The preceding results indicate that linearized theory fails 
to describe properly the propagation of sound in the vicinity of 
a near-sonic throat. As has been discussed previously (refs. 2-5), 
the size of the region around x=0 in which the outer expansions 
fail is of 0(e). Therefore a "stretched" or inner variable 


f (y+1) a . \ x 
1 2 J e 


(4.1) 


3 ^ 

and a modified time variable T = a 2 t are introduced, and the solu- 
tion to the systems (3.2), (3.4), apd (3.6) is examined in the 

inner region, X=0(1). The factor a^ in X and T is included for 
convenience in order that a does not appear explicitly in the sub- 
sequent analysis. The basic flow Mach number in the inner region 
has the expansion (refs. 1,5) 


M ( x ; e ) = M 1 (X ; e ) = 1 - en^ (X) + 0(e 2 ) (4.2) 

2 ^ 

in which m-^ (X) = (1 + X ) 2 . The perturbation quantities are ex- 
panded in the form 


p(x,t;e) = p 1 ( X , T ; e ) = er 1 (X,T) + e 2 r 2 (X,T) + ..., 

u (x, t; e ) = u 1 (X,T;e) = ey 1 (X,T) + e 2 y 2 (X,T) + ..., (4.3) 

p(x,t;e) = p 1 (X,T;e) = ea 1 (X,T) + e 2 a 2 (X,T) + ... 


Substitution of the expansions (4.2) and (4.3) into the system 
(3.4), changing the spatial variable to X, and equating the co- 
efficients of each power of e to zero following the same steps as 
in reference 2 yields at the lowest order 


% 

+ y-, =o, 

X 

(4.4) 

(y x 

- i Oj) + (2p x + r x - - 2m 1 )p 

J- X 

(4.5) 


2m' 

"y+1 (2lJ l + r l " a l } = °' 


a i 

X x 

- yr =0. 

X 

(4.6) 
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For later purposes it will be necessary to use the second order 
counterpart of the inner equation of continuity (4.4), which is 


Mr 


= -r. 


" y l r l 


x 




(4.7) 


The corresponding inner form of the shock conditions (3.6) 
can be obtained using the expansions (4.2) and (4.3) and by rep- 
resenting the unknown shock location in the form 


T = H(x;e) = H Q (X) + (X) + ... . 


Since t and x are scaled according to L/c s and L, 
this yields the expression for the shock velocity 
be determined in the course of the analysis, as 


V = ( 


2 ec s 

Y+l Hq ( x) 


+ 0(e ) . 


respectively , 

V, which is to 

(4.8) 


Equations (4.2), (4.3), and (4.8) are substituted into the equa- 

tions (3.6) and the coefficients of each power of e equated to 
zero to yield the lowest order shock conditions on the inner solu- 
tion in the form 


[r l + y l ] = °' 

[r^ 2y^ y ^1^ = 

tv i - 7 + “o'T 1 “i 


r i U l + ° 1 V 1 + 2m i y i ] 


0 . 


(4.9) 

(4.10) 

(4.11) 


The quasi-linear system (4. 4-4. 6) describes the behavior of 
sound propagating in the near-sonic duct throat region to leading 
order in the small parameter e. It expresses the fact that the 
sound propagation process is inherently nonlinear in the near- 
sonic flow. Solutions to this system will be discussed in the 
following, where it will be seen that in general discontinuities 
in the quantities r^, y^, and cr^ will develop. The magnitude of 

these discontinuities as well as their unknown location T=Hq (X) 
will be determined using the shock conditions (4.9-4.11). 


Boundary Conditions 

The theory outlined in the preceding sections will be applied 

in the following to study the sound field propagating upstream in 

a Crocco-Tsien duct. For convenience, the upstream exit of the 

duct at x = x <0 will be assumed anechoic. It is known (ref. 3 ) 
e 
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that vanishing of the downstream propagating wave at a point in 
the outer region requires M Q u 0n + p Qn = 0 there. Thus, the ane- 

choic end condition will relate the complex constants A^ n and 

B in the outer solution (3.11) in x<0 according to 


B ln 

A In- 


V x e )G 3n (x e )+G ln (x e> 

M n (x^) G 4n (x^) +G 0 _ (x,J 


2n 


(4.12) 


The superscript on the constants in equation (4.12) is used to 
indicate that these are the constants appropriate to the linear 
solution (3.11) in x<0. On the downstream side of the throat the 
constants will, in general, be different and will be denoted 
A-^, when necessary. 

The remaining boundary condition is an acoustic source con- 
dition. Here it will be assumed that the source is located at 
x=Xq, and that the perturbation velocity there is specified, i.e., 

U(x;e)u(x Q ,t;e) = c Q (x Q )M 0 (x Q ; e) u (x Q ,t ; e) = S sin tot 

where S/cq(xq)<<1. Two cases arise according as the point x=Xg 
lies in the inner region, x 0 =0(e), or in the outer region, 

Xq=0(1). In the former case, the source amplitude is required to 
be proportional to e , so that the boundary condition becomes 

c 0 (eX 0 )M 1 (X 0 ;e)u 1 (X 0 ,T;e) = eU sin QT, 

^here the modified frequency fl=to/a 2 has been ^Introduced and where 
U is a constant. The inner point Xp=-T (y+1) a ^XQ _ p(l). Using 

Equations (3.5), (4.2), and (4.3), the source condition on the 

leading order inner problem becomes 

^(X^T) = -Ij- N sin QT (4.13) 

v^here 

N = (X±i) 3/2 _U_ 
y 2 ' C 

s 

On the other hand, if the acoustic source is located down- 
stream of the throat in the outer region the amplitude S will be 
taken to be of 0(6), consistent with the definition of 6 as the 
perturbation level in this region. In this case the source 
condition becomes 


16 


c 0 (x 0 )M(x 0 ;e)u(x 0 ,t;e) 


6c 0 (x 0 )M 0 (x 0 )u 0 (x 0/ t) + . . . 


6U sin tot. 


In this situation the boundary condition on the leading order solu- 
tion in the outer region is 

u 0 (x 0 ,t) = N sin tot (4.14) 


in which N = U/cq (xq)Mq (xq) . Using the solution (3.11) in x>0 , 

the condition (4.14) yields relationships between the complex 
constants Aj+ n and B-j^ in the form. 


B 


+ 

In 


- R+(x )A 
n o' 


+ 

In 


F 

n 


f 


where the complex constants F n =0 if n^l, and is proportional 
to N. It is easily verified using Equation (3.11) that 

R n (x 0> = G 3n (x 0 ) / G 4n (x 0 ) * 

In the specific results derived later the major emphasis is 
on the problem for which the acoustic source lies in the inner 
region, Xg=0(e), so that the source condition is as in Equation 

(4.13) . The reason for this emphasis is that the throat source 
problem can be solved analytically. Thus it offers an overwhelm- 
ing advantage for gaining insight into the nonlinear effects of 
the near-sonic flow on sound. In reference 4 the authors discussed 
the more practically interesting latter case, for which condition 

(4.14) is relevant, under the restriction that shocks do not form. 
The complete solution in this case involves a complicated iteration 
scheme. The generalization of this scheme to include the develop- 
ment of shocks will be the subject of a future publication. 
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5. SOLUTION OF THE INNER PROBLEM: 
ACOUSTIC SOURCE NEAR THE THROAT 


The solution of the system (4. 4-4. 6) subject to the shock 
conditions (4.9-4.11), the boundary condition (4.13), and appro- 
priate asymptotic matching conditions to be discussed later can be 
effected most easily by introducing certain physical considerations 
to reduce the complexity of the problem. 


Reduction to a Signaling Problem 

It has been shown in references 3 and 5 that the combination 
r l + Pi is the leading order inner representation of the downstream 
propagating Riemann variable for the system (3.4). It satisfies 
the linear relation (r^ + y-]_) x = 0 and is continuous in the inner 
region by virtue of condition (4.9). This, together with the fact 
that the combination Pq + uq of the outer solution remains finite 
as x-^-0 (see eqs. (3.11) and (3.12)), suggests that the downstream 
propagating wave is unaffected at the lowest order by the nonlin- 
earity in the inner region and that, in fact, r^ + y^ vanishes. 

Also, since discontinuities in the inner solution represent 
shocks whose strength is 0(e), it is anticipated that entropy 
changes occurring near the throat are of 0 (e-^) . If so, the first 
order inner perturbation field is homentropic, and the appropriate 
solution of equation (4.6) is = yr]_. 

Thus, it will be assumed here, and later verified by asymp- 
totic matching, that y]_ =-ri and ai = yr^. In this event, equa- 
tions (4.4), (4.6), (4.9), and (4.10) are satisfied identically. 

The remaining inner equation of motion (4.5) and shock condition 
(4.11) reduce to 


2y + { (y+1) y, - 2m. }y. - 2m' u. = 0, 

1 T X 


[2y. ] + HA [- 


(y+1) 2 


y 1 + 2m 1 y 1 ] = 0, 


Finally if the definition 


(X,T) = ^ y x (X,T) - m 1 (X) 


(5.1) 


(5.2) 


is introduced, the inner problem (5.1) , (5.2) , and (4.13) is trans- 

formed into the following signaling problem: 


n T + nn x = x, -» < x < x Q , 


(5.3) 


n(X 0 ,T) = -(1 + Xq)' 2 + N sin fiT, 


(5.4) 
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where discontinuities in II occur on the space-time curves T=Hg(X) 
and satisfy 

2 [n] - Hg[JI 2 ] = 0 (5.5) 

Interest here is directed toward time-periodic or steady-state 
solutions. Hence, no initial conditions are stated, i.e., it is 
assumed that the source has been acting for a long period of time. 

Equation (5.3) can be solved using the method of characteris- 
tics (ref. 11) it states that 


dn 

dT 


(5.6) 


along the characteristic curves in (X,T) which satisfy 


dX 

dT 


n 


(5.7) 


A parameter x is introduced such that x = constant along a charac- 
teristic curve, and x is chosen such that x=T when X=Xq . Differ- 
entiation of equation (5.6) with respect to T holding x fixed and 
use of equation (5.7) leads to 


dT 

along a characteristic. Hence II is of the form 

X = A(x) sinh T + B(x) cosh T. 

Application of the boundary condition (5.4) and the condition T=x 
on X=Xq then yields the inner solution in the form 

IT = <}>(x) cosh (T-x) + Xq sinh (T-x), (5.8) 

X = <j>(x) sinh (T-x) + Xg cosh (T-x) , (5.9) 

in which 

4>(t) = N sin fix - (l+Xg)* 2 . (5.10) 

Equations (5. 8-5. 9) are a parametric representation of the 
exact solution to the inner problem (5.3) , (5.4) . As will be seen 

in the following, this representation can be matched asymptotically 
to the outer solution even though it is not possible to express 
II explicitly as a function of X and T. 

The behavior of the inner solution (5.8) and (5.9) when Xg=0 
has been discussed extensively in references 3 and 5 under the 
restriction that II remains continuous over -°° < X <0. Here, 
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however, solutions which develop shocks are of primary interest. 
Shocks occur in the solution n whenever the one-parameter family 
of characteristic curves (5.9) overlap in the (X,T) plane, thereby 
rendering the expression (5.8) for II multivalued. The condition 
for the occurrence of such multivaluedness follows from the fact 
that the boundary of a region of overlap must be an envelope of the 
family of characteristic curves. An envelope, if one exists, is 
given by simultaneous solution of equation (5.9) with its deriva- 
tive with respect to x : 

<Mt) cosh (T-x) - (<J>'(t) - X Q ) sinh (T-x) = 0. (5.11) 

Now, in the present problem, only those characteristics pro- 
pagating upstream from X=Xg can carry the solution through the 
throat region; i.e. dX/dT < 0 at X=Xo on all characteristics of in- 
terest. Equation (5.7) then indicates that n<0 on each of these. 
Since n=cf> at X=X 0 it follows that cj)<0 on every characteristic prop- 
agating upstream from the source. Now, T-x>0 on any such charac- 
teristic, and hence equation (5.11) implies that cfi ' — X q < 0 and that 
<f>/ (<f> '-Xg) <1 on all characteristic curves which overlap and form an 
envelope. The second of these conditions provides the criterion 
for the development of shocks in the inner solution II ; shocks will 
form if and only if 

(x) - X Q < cj) ( x ) (5.12) 

for some range of x between, say, x = 0 and x = 2TT/f2. The inequality 
(5.12) is the condition on source strength, frequency, and location 
such that the inner solution II will develop shocks. 

Using equation (5.10), the condition (5.12) implies 

N ( sin fix - fl cos fix) > (1+Xq)* 5 - X Q . 

If the left side of this inequality is maximized with respect to x 
it follows that it can be satisfied for some range of x over each 
period if and only if 

Nd+fl 2 )^ > (l+X 2 )* 2 - X Q (5.13) 

Equation (5.13) is the specific relation for the present problem 
which determines in terms of the given parameters whether the non- 
linear inner perturbation field develops discontinuities. It was 
derived for the special case Xq=0 in references 2,5 and 6. 

Wherever N, fl, and Xq satisfy condition (5.13), the upstream 
propagating wave represented by n will steepen and break in typi- 
cal nonlinear fashion. Breaking first occurs at the point (X,T) 
where the envelope of characteristics forms a cusp. The charac- 
teristic x through the cusp is determined by eliminating T between 
equation (5.11) and its derivative with respect to x; 

(2<j) 1 (x ) -X q) cosh (T-x) - (cj>"(x) + ( x ) ) sinh (T-x) = 0. 

Consideration of the simultaneous solution of this equation along 
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with equations (5.9) and (5.11) shows that, for N slightly greater 
than the critical value indicated in equation (5.13), the cusp is 
located at large negative X. Also, using the same three equations 
it can be seen that the cusp will lie on X=X n only if <f) ( t ) = 0 and 
2(f>' (t)-X 0 =0, i.e., if u 

2 J- 

N sin fix = (1 + Xq) 2 and 2Nfi cos fix = Xg , 
which requires that 

N = N c = (1 + Xg + Xg /4a 2 ) 1 *. (5.14) 

With the information developed thus far it is possible to 
establish clearly a qualitative picture of the development of 
shocks in the perturbation field arising from nonlinear interac- 
tions near the duct throat. In general, the inequality (5.13) 
indicates that shocks will form as any of the parameters N, Q, and 
Xg increases. For the present discussion ^ and Xg will be consid- 
ered fixed. Then for sufficiently small N the solution is contin- 
uous, as in references 3 and 5. As N is increased beyond the value 
indicated in equation (5.13), a shock wave occurs in each period 
of the inner solution IT. Provided N is less than the value N of 
equation (5.14) the shock starts with zero strength at a finite 
distance upstream of X=Xg, where the steepening wave profile first 
develops an infinite slope. As N increases towards N c the shock 
formation point moves back towards X=Xg. For N>N C no cusp forms 
on the envelope of characteristics in X<Xg. In this case, as will 
be discussed in more detail later, the shock forms directly at the 
source, X=Xg. 

The actual characteristics over 0<x <_2it /Q , given by equations 
(5.9) and (5.10) for the present problem, are shown plotted in 
figures 2(a) to 2(d) for the case Xg=l, fi=l. In this instance, 
the shock exists for N>.222, and the value of N c is 1.5. These 
figures will be discussed in detail later, but for the present it 
is of interest to note that as N increases towards N c the cusp 
turns so that the slope of the envelope (and thus of the shock) 
becomes positive there. That is, the shock actually travels down- 
stream towards the source after its inception. Because of this, 
there is a range of N slightly below N c in which the shock forms 
with zero strength at the cusp, travels downstream to interact 
with the source, and than proceeds back upstream through the 
throat . 


Shock Fitting 

The multivaluedness predicted by equations (5.8) and (5.9) is 
resolved by the introduction of appropriate discontinuities which 
satisfy the jump relation (5.5). The approach to be used in this 
process is parallel to that of Whitham (ref. 10) . However, the 
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inhomogeneous nature of the differential equation (5.3) signifi- 
cantly complicates the analysis, owing mainly to the fact that the 
characteristics are not straight and the solution II does not re- 
main constant along them. Fortunately, however, it is possible to 
carry out the shock analysis analytically despite the complications. 

Each point of a shock path in space-time is a point at which 
two characteristic curves, carrying different values of n, inter- 
sect. Let the characteristic intersecting the shock from ahead 
(i.e. from larger negative X) be labeled T? and that intersecting 
the shock from behind be Xj_. Then x^>X 2 on the shock. Using 
equation (5.9) , a point X on the two characteristics, and on the 
shock, satisfies 

X = Xq cosh (T-x^) + <f>^ sinh (T-x^) , (5.15) 

X = Xq cosh (T-x 2 ) + <f >2 sinh (T-x 2 ) , 

in which <j>^ = ^(x^), <j> 2 = 4> ( x 2 ) . The shock condition (5.5) gives 


dT 

dX 


Hq (x) 


n i +TI 2 


(5.16) 


where = n(T;xj.), II 2 = n(T;x 2 ). The three equations (5.15-5.16) 
determine simultaneously the parameters x^(X) and X 2 (X) and the 
path T=Hq (X) along the shock. 

It is now required to derive a generalization of the equal 
area relation of weak shock theory using equations (5.15-5.16). 
First, equations (5.15) are differentiated with respect to X along 
the shock and added, which leads to 

dT 

2 = ^ {Xq sinh (T-x^) + <j>^ cosh (T-x^) 

Xq sinh (T-x 2 ) + <f >2 cosh (T-x 2 ) } 
dx 

+ (~X 0 sinh (T-x x ) + 4*^ sinh (T-x^ 

-<j>, cosh (T-x^) } 

dT 2 

+ {-X Q sinh (T-x 2 ) + §' 2 sinh (T-x 2 ) 

— 4> 2 cosh (T-x 2 ) } . 
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The term multiplying 


equation (5.8) to 
then reduces the result 


dT 
dX 
be IIj 


in the above 

+ II?; use of 
to 


relation is seen using 
the shock condition (5.16) 


dx. 


dx 


dx 


dx . 


n 


1 dX 


+ n 2 = ((>1 sinh (T-t 1 ) + $' 2 sinh (T-x 2 ) (5.17) 


It is to be shown that equation (5.17) can be integrated to 
yield an equation for x^(X) and X 2 (X). To do this requires first 
that the unknown function T(X) be eliminated. This can be accom- 
plished by using equations (5.8) and (5.9) to get 


sinh (T-x) = 


cosh (T-x) = 


xd)-x 0 n 

,2 2 

<f> x Q 

nct>-xx 0 

,22 
4> -x Q 


(5.18) 


(5.19) 


Equation (5.18) is now substituted into equation (5.17) , and the 
equation is multiplied by n 1 - n 2 (which is non-zero in general) , 
to obtain 

2 dx ;i 2 2 

<n l- n l n 2 ) dT + ("iW dx- 


dx 


= x(n^n 2 ) { 


A 2 V 2 dx 

*i‘ x o 


1 , *2*2 ’ dT 2 , 
' **% XI TT J 


i 2 v 2 dX 

V X 0 


2 1)5 1 

- X 0 (n i“ n i n 2 ) { T i 


dr. 


^2 


dx , 


v 2 dx 

V x 0 


2 2 

^2- x 5 


dX 


Using equations (5.8) and (5.9), it can be seen that 

2 2 2 2 
IT - YT = <j> - Xq . 


(5.20) 


(5.21) 


Thus the left side of equation (5.20) can be written as 

99 9 dX 1 9 9 9 dx 

{ *1 " x q " (II i n 2 " X )} ~3T + {(II i n 2 " X > ~ *2 + X 0 } 


dX 


while the second term on the right side is 
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dx 


-x«<K 


0 Y 1 dX 


1 , X 0 (n 1 II 2- x2 > 

+ — <K 


dx . 


2 2 
^l- X 0 


1 dx 


dx. 


+ x rt 4>‘ 


0 r 2 dX 


x Q (n 1 n 2 -x ) 

.2 v 2 

^2 X 0 


dx. 




2 dX 


Use of these two results allows equation (5.20) to be rearranged 
in the form 


j. 2 dT l j. 2 dT l _ v 2, dT l 

9-^ “* rn ** r\ \ 


2 dX 


O' dX 


dx„ dcf>, 

-)-x n (- 


d<J). 


dX 


k 0 ' dX 


dX 


4>t 4>i dx, 4> 9 4>A 

+ X(n-n 2 ) { 2 2 -gx" + 2 y 2 -dx" } 

♦l " X 0 V X 0 


2 dT 1 

+ (n n 2 -x H-jjj- 


dT 2 + Vi dT 


dX 


.2 v 2 dX 
♦l _X 0 


( 5.22) 


Vi dT 2 -l 

.2 „2 dX 

V X 0 


that 


Now, it is easily shown, using equations (5.18) and (5.19) 
cosh (x 1 -x 2 ) = cosh {T-x 2 - (T-x^ } 

(n 1 n 2 -x 2 ) (<t) 1 cj) 2 -x 2 )-xx 0 (<}, 1 -(j) 2 ) (n ] _-n 2 ) 

” 9 2 2 2 ' 

^ l - X 0 } ( * 2 "V 


and, similarly. 


sinh (x^-x 2 ) = 


x(n 1 II 2 ) ( 4 s 2_ ^ 2 V x q ( 4 *i 4 * 2 ^ X ^ 

^2- X 0> 


These two relations are considered as simultaneous algebraic equa- 
tions for the quantities Xd^-I^) and (n i II 2 -X 2 ) whose solution is 


X(n 1 -n 2 ) = (4> 1 4> 2 -Xq) sinh (x 1 -x 2 ) + X Q (<f> 1 -4> 2 ) cosh (x^x^ 
n i n 2 _x2 = ( 4* 1 4> 2 — X 0 ^ COSh (t 1 _X 2 ) + X 0 ((f> l" <f) 2 ) sinh (t 1 _X 2 ) 


(5.23) 
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Finally, equations (5.23) are substituted into equation 
(5.22). When this is done it is a matter of straightforward alge- 
bra to show that the coefficient of sinh (x^-x 2 ) simplifies to 

dx.. dx„ dx. dx 9 

V 2 "dX* + ^1^2 ~3X + X 0 ((f) l“ <,> 2 ) ( "dX diT 5 ' 

while the coefficient of cosh 




, dT l 


^0 T 1 dX 0 T 2 dX 

Equation (5.22) then becomes 


dx~ „ dx. dx~ 

x n<f>o -3v- + ( < f ) i < f ) 2 “ X 0* (_ dX 'dX _) ’ 


^1 dX ^2 dX X 0 dX ^ T 1 T 2 ^ Xn 


k 0 dX ""1 T 2 ; 


+ ^dX ^1^2^ + X 0 dX ( t 1 -t 2^ sinh ^ T i x 2^ 

+ {X 0 dX ^i“^2 ^ + ^1^2 “ X 0> dX ( T 1~ T 2^ cosh ^ T l _T 2^ 
which can now be seen directly to be 


dX / * dT X 0 dX (T 1 T 2 ) X 0 dX ((J> 1 <t> 2 ) 


j 2 

4" { (^^^2 x o ) ( t^-T 2 ) ^ cosh 


Hence, equation (5.17) is equivalent to 
x. 


f <p 2 dx = Xq (x 1 ~x 2 ) - Xq ( 4> — ^ 2 ) + ^1^2 “ X 0^ sinh (t^-t 2 ) 

T 2 

+ X Q ( 4> 1 — cf> 2 ) cosh (x ^-x 2 ) + D, (5.24) 
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where D is a constant of integration to be determined presently. 
Equation (5.24) is the generalization to the present problem of 
the so-called equal area relation of weak shock theory (ref. 9 ) . 

It is a transcendental relation between Ti (X) and x 2 (X) on the 
shock, which, together with equations (5.15), completely determines 
the location and strength of the shock in the inner region. 


Analysis of Characteristics and Shock Formation 

The determination of the constant D (eq.(5.24)), which depends 
on the parameters N, ft, and Xq, requires a detailed analysis of 
he configuration of the characteristic curves. These are illus- 
trated (over one period) in figure 2 a-d for increasing values of 
N with XQ=ft=l. The circular symbols denote points on the envelope, 
calculated using equations (5.9) and (5.11). It can be seen from 
equation (5.9) , and also from the figure, that so long as 
-4 >(t)<Xq, there exists a set of characteristics whose slope de- 
creases towards negative infinity at some X<Xg at a time which is 
easily found from equation (5.9) to be T=x +tanh - !(-<(> (x ) /Xq) . Each 
of these characteristics is a hyperbolic cosine curve, symmetric 
about its turning point, which therefore returns to intersect the 
source line X=Xq at 


T = x + 2 tanh 1 (-<(> (x ) /X Q ) . (5.55) 

Equation (5.25) can also be found directly from equation (5.9) by 
setting X=Xq and using the identities relating sinh z and cosh z 
to tanh hz . The value of II, which vanishes at the turning points, 
becomes positive on the downstream propagating portion of the cosh 
curves . 

It can be shown that the combined Mach number U(1 +u)/cq is 
equal to 1+eII in the inner approximation. Hence, when a charac- 
teristic assumes a positive slope it means that the net flow 
(steady flow plus acoustic perturbation) is supersonic. The turn- 
ing of the characteristics seen in figure 2 therefore represents 
an actual acoustic choking of the duct flow whereby sound imparted 
to the flow during portions of the period when the source amplitude 
is sufficiently low is prevented from propagating upstream through 
the throat. This choking phenomenon significantly complicates the 
analysis of the development of shock waves in the duct. 

For N sufficiently small, the shock must first form at the 
point in X<Xq at which the envelope of characteristics forms a 
cusp (see figures 2 (a) and (b)). At such a point, x,=T 2 must be a 
root of equation (5.24), which implies that D=0. As seen in figure 
2(b), however , when N increases towards N c the shock assumes a posi- 
tive slope at its inception and initially travels downstream 
towards the source. Provided that it turns back upstream before 
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reaching X=Xq the value of D remains zero. The condition that 
the shock turn exactly at X=Xq requires that X]_=T, and that 
n(T 1 )+n(x 2 )=0. It follows from equations (5.25) and (5.8) that 
n(T 2 )=-4>(T 2 ) at X=Xq , so that <f> (x^) =cf> (t 2 ) . Thus 


x£ = Tr/ft - x* 

x* = x* + 2 tanh -1 (-<j> (x*) /X Q ) 

are the two equations which determine the characteristics x|, x 2 
intersecting on the shock when it reaches X=Xq with zero velocity. 
But these values also satisfy equation (5.24) with D=0, which 
yields the upper limiting value of N=N* for which D=0 . It does 
not appear possible to derive a simple analytical result express- 
ing N* in terms of Xq and ft. However, numerical solution of the 
two above equations together with equation (5.24) in the case 
Xo=ft=l yields the result N*=l. 12717 (the value of N c in this case 
is 1.50). The general configuration of the shock near X=Xq is 
indicated schematically in figure 3 (a) and (b) , where the envelope 
of characteristics is shown by dashed lines, and the shock is in- 
dicated by a heavy line. The picture, of course, is repeated pe- 
riodically in T. 

In the following discussion the form assumed by equation 

(5.24) when the shock intersects the source line X=Xg will be 
needed. At any time at which the shock is at Xq the value of x-^ 
equals T and the value of x 2 is that appropriate to the returning 
branch of the cosh characteristic at that T; i.e., x„ (see eq. 

(5.25) ) is given by 

x n = T 0 + 2 tanh (-4> 9 /X n ) . (5.26) 

cf >2 and X=Xq yields 

(5.27) 


(5.28) 


If equations (5.27) and (5.28) are substituted into equation (5.24) 
it becomes 


Use of equations (5.18) and (5.19) with JI^^ - 


sinh (x^-x 2 ) = 


cosh (x^-X 2 ) = 


2 v 2 
* 2- X 0 ' 

* 2- X 0 
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(5.29) 


T 2 

/ ( 4> 2 — Xq ) dr = 2X 0 (j> 2 + D. 

T 1 

Equation (5.29) is the form of the generalized equal area relation 
(5.24) when the shock is at X=Xq. 

2 J- 

For N* < N < (1 +Xq) 2 the upper branch of the envelope moves 
towards X=Xg until u it actually touches the source line when 
N=(l+Xy) 3 s. The details of the characteristic picture are not vis- 
ible on the scale of the actual plots used in figure 2 . They are 
shown schematically in highly exaggerated form for the case 
N=(l+X§) J s in figure 3(c). Over this range of N the shock starts at 
the cusp, travels downstream, and intersects the source line with 
a nonzero velocity. Thus, the shock splits into two branches. 

On the lower branch, which includes the cusp, the value of D in 
equation (5.24) is zero. The characteristics t + and x_ indicated 
on the figure are the two which intersect on the lower branch of 
the shock at X=Xq . They are determined by simultaneous solution 
of equation (5.28) with D=0 and equation (5.26) using t^= x + and 
T 2 =t_. It is noted that this lower branch of the shock is actually 
very short. In the case illustrated in figure 2, Xq= 1 and the 
cusp is located at about X=0.999. 

Directly above T=t + the solution n is determined on X=Xg not 
by the source data but rather by the values appropriate to the 
returning cosh characteristics x£x_. Over this region the trans- 
formation relating t to T is given by equation (5.25) . Here the 
acoustic choking of the flow has extended all the way back to the 
source line X=Xg. The source, whose amplitude is insufficient to 
overcome the oncoming supersonic flow which it experiences, is 
completely blocked from imparting sound in the upstream direction. 
This suggests the physical argument which determines the value of 
D appropriate to the upper branch of the shock, which is the branch 
which actually travels back upstream towards X=-°°. The shock must 
begin at the point T=t£ at which the source amplitude <J> (which is 
negative) exactly cancels the value of II on the returning charac- 
teristic x£ to yield sonic conditions in the flow and a shock 
velocity of zero. This occurs just as in the previous determina- 
tion of N* only if <J> (x£) =<j> (x|) and if equation (5.26) is satisfied 
by x£ and x^. Thus, 


x£ = it /Q. - 


x 


2 ' 


T° = x°2 + 2 tanh 1 (-cf)°/X 0 ), 
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where <J>°Ecj) ( t °) , are the two equations which determine the charac- 
teristic curves intersecting at the initial point of the upper 
branch of the shock. The physical argument above is supported by 
detailed numerical studies performed by the authors which show 
that there is one and only one such value t£>7t/2S2 at which the 
above two equations are satisfied. 

O O 

Once t 2 . and X2 are known from the above equations then D is 
determined for the upper branch of the shock by using equation 
(5.29): 

O 

T 1 

D = / ( <j) 2 — Xq ) dx - 2X Q (f>2 • (5.30) 


For the range (1+Xg) 3s <N<N c the characteristics appear as in 
figure 2(c). A schematic representation of their behavior near 
X=Xg is shown on figure 3(d). The only difference between this and 
figure 3(c)is that there exists a range Xg<T<x^ over which <{>>0, 
where 

, 1-Xg + 

x“ = sin" ( -fr ■) 2 and Xg = ir/ft - x Q . 

The envelope of characteristics has split (see figure 2(c)), but 
there still exists a cusp as indicated in figure 3(d). The values 
x + , x_, x£ and x£ are determined exactly as above, and D is given 
for the upper branch of the shock by equation (5.30) . 

For N>N C the cusp no longer exists, as was described earlier; 
the lower branch of the envelope, and hence the shock, disappears. 
This is illustrated by the plot of figure 2(d). Although no corre- 
sponding schematic is shown in figure 3, the upper branch of the 
shock starts at X=Xq exactly as illustrated in figure 3(d). The 
values of x£, x£ and D are determined as above. 

In the asymptotic matching to be discussed in the following 
chapter, the value of D on the shock which propagates to large 
negative X is required in order to use equation (5.24) to find 
the limiting values of x^ and x? as X-*- °°. This is what has been 
determined in the preceding rather lengthy analysis. For the 
remainder of this work the integration constant for equation (5.24) 
will be denoted D (N) , where D(N)=0 for N<N* , and D (N) is given by 
equation (5.30) for N>N*. 

For the special case Xg=0 it is easily shown that N*=l and 
that equation (5.30) yields 0 


29 



ftD(N) 


= (1+ ^-) (Tr-2sin 1 ^•)-3(N 2 -l) 35 

for N>1 . In this case no lower shock branch ever appears. 

This result was derived previously in reference 6. The numerical 
examples to be discussed later are restricted to this special case. 
The shock relations are used in the form of equation (5.24) and 
the first of equations (5.23). For a fixed X the two equations 
are solved using a Newton-Raphson procedure. Once x-^(X) and 
t 9 (X) are known then T=H n (X) is calculated using one of equations 
(5.15). 0 
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6. ASYMPTOTIC MATCHING 


The solution to the inner problem has been completely deter- 
mined in the preceding chapter. This solution must be asymptoti- 
cally matched to the outer solution to determine the perturbation 
level P of the acoustic field which propagates upstream from the 
throat and the constants A^ and B ln of equation (3.11) . The 
matching will be carried out on the perturbation velocity u; cor- 
responding results for P and p follow. 

The one-term outer expansion of u is u- L (x,t;e). A one-term 
inner expansion of this follows by introducing the inner variable 
X so that u ~ 6u^(x(X;e) , t;e) , and expanding to leading order in 
e keeping X fixed. This is equivalent to a leading order expan- 
sion of u^ for small x and e, which is precisely the limit uq 
given by equations (3.11) and (3.12). Thus 

u ~ -| — r E A exp {into (t - — p- £n |x|)}. 

I x l n=0 xn a 2 

Finally, the physical value of u is the real part of this expres- 
sion. Writing Aj = A n exp (intoi^ n ) , this yields for the one-term 
inner expansion of the one-term outer expansion of u the expression 

u ~ -I — p E A„ cos nto ( t — — p" &n | x | + 4^ ) • (6.1) 

|X1 n=0 a 2 

The one-term inner expansion of u is, from equation (4.3), 


u 


ey 1 = {(1 + X 2 ) h + n(X,T)} 


A one term outer expansion of this is derived by using equations 
(4.1) and (5.21) and expanding each term to leading order in e 
keeping x fixed. Thus 2 2 

n r . 2 „2 2 h r Y+l _ n 3s Ix| n * " X 0 e 2 + . 

n - - {( J) - X Q + X } -{— a > £ u + (y+l) a ^2 


} 


(l + x 2 ) 1 ' 2 


which yields 


{ I+1 ^ M (1 + 


u 


O L4. 

" x | 2 


rY+l 


a } 


-3/2 


(y+l) a x 2 


+ ... } , 


{l - (<j> 2 - X 2 ) } 


( 6 . 2 ) 
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In conjunction with this result, the outer representation of 
the characteristic curves can be determined by similarly expanding 
equation (5.18). This results in 

T = a^t = x-£n|<j)(T)+X 0 |+g(x;e) (6.3) 

to leading order, where 


g ( x ; e ) 


£n 2 ( 


Y+l 

2 




It is noted that in deriving equation (6.3) use has been made of 
the fact that as X->— 00 both II and <j> are negative on the character- 
istic curve. Since (T-t)>0, equation (5.18) implies that <J> 2 -Xg>0. 
It follows that 4 >+Xq< 0 on each characteristic which extends to 
large negative X. 

Equation (6.2) is the one-term outer expansion of the one-term 
inner expansion of u expressed in parametric form, where the param- 
eter t is defined in terms of x,t by equation (6.3). 


Now, because cj> ( x ) is a periodic function of t (see equation 
(5.10)), and T-t is a periodic function of t, it follows that <}> 2 (t) 
is a periodic function of T with period 2ir/Q and thus of t with 
period 2 tt/o). Hence <j> 2 * (x) can be written in the outer region as a 
Fourier series. It is convenient to include Xq in the expansion 
and write 

2 2 C 0 (x) 

d) (t) - X_ = — ~ + £ C (x) cos nwt + D (x) sin nut 

Y 0 2 . n n 

n=± 


where 


2ir/o) 


jc \ bi r . 2. > 2. /cos ncot\ 

= ? J <t) - V( sin mot) 


dt 


2n/Q 


= - f ( cj) 2 ( T ) - X 2 ) | C ? S Jox} dT 

7T J vy ' ' 0 | sin nflT/ 


(6.4) 


To carry out an explicit evaluation of the Fourier coeffi- 
cients it is necessary to consider the shock relations (5.23) and 
(5.24). In the outer limit, as X-*- ■ «>, it is easy to show from 
equation (5.21) that 


x(n 1 -n 2 ) 


2 2 
♦ l "*2 
2 ' 
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so that the first of equations (5.23) becomes 
2 2 

* l -*2 2 

2 = (4> 1 4>2 “ x o ) sinh (t i” t 2^ + X 0 (<fl l _<() 2 ) cosh (T 1 -T 2 ) (6 - 5) 

It follows that equation (5.24) takes the form, as X->- <», 

2 2 

2 4>I"*2 

X 0 (x 1" T 2 ) “ X 0 (<J> 1 -<|) 2 ) + 2 ~ + D(N) • (6 * 6) 



The two transcendental equations ((6.5) and (6.6)), which are in- 
dependent of X, determine the limiting form of t^(X) and T 2 (X) on 
the shock wave as it travels out of the throat region. These two 
equations were solved numerically using a Newton-Raphson scheme 
in computing the results to be presented later. 


Assuming that the outer limits of x^ and x? are known from 
solving equations (6.5) and (6.6), it is possible to proceed with 
the calculation of the Fourier coefficients C n and D n . The proce- 
dure to be used is similar to that of Blackstock (ref. 12). Since 
<J> is known explicitly as a function of x , it is desirable to change 
the variable of integration in equation (6.4) to x . In doing so, 
however, it must be borne in mind that the mapping from T to x 
becomes discontinuous in the presence of shocks in the sense that 
over any period T 2 <T<T 2 + 2tt/£ 2 at a given X, the values of x range 
only over x^<x<X 2 + 2ir/fi, where X 2 = x (X,T 2 ) and X]_ = x(X,Tj). The 
characteristics T2 <t<t 1 can ke re 9 ar cled as having been "eaten up" 
by the shock. Thus, after shifting the integration range to 
T 2 <T<T 2 + 2t\/Q, equation (6.4) can be written 



(4> 2 (x) 


y 2. i cos n^T(x) ) dT 
A (r (sin nf2T(x) | dx 


in which T(x) 
<j>2 - Xq as ($ 
the form 


is as given by equation (6.3). Writing the factor 
+ X 0 ) 2 - 2 Xq ( $ + Xq) , these integrals can be put in 
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n 


c + iD = 
n n it 


x 2 +27r/fi 

J (*tx 0 ) 2 a in8T(,) |dt 


-2X, 


§ / 


t +2ir/fi 
(• z 


,, , v v infiT(x) dT , 

(<p + X Q ) e -pjr dT. 


(6.7) 


The first integral in equation (6.7) becomes, after using 
equation (6.3) , 


/ 


T 2 + 2TT/f2 


<* + V 2 e inOT(T) (1 ' dT - 


( 6 . 8 ) 


which, upon integrating the term containing ' by parts, can be 
written 

x 2 +2ir/fi 

„ >2 infiT(x) , , 1,, , „ ,2 infiT(x..) 

( 4> + X Q ) e dx + 2^\ + X o J e 1 


1 flk . v \ 2 infiT ( 

" 2 (<t> 2 + V e 


T2 + 2n/tt 


x 2 ) + insy* (<() + Xq 


,2 infiT (x) dT , 

) e a? dx. 


The last term of this is seen to be proportional to the original 
integral (6.8); thus that integral can be written, after substi- 
tuting equation (6.3) for T(x), as 


1+infi 

2 infig T 

2~2 e 1 

HnV 

4 


In' 


where 

T ln 


x 2 +2 tt/ fi 


v infix 


= J (<(,+X 0 ) v e J - li “ l dx + 2 


ki U 1+ xJ V e in ^ T l-U 0+ xJ V e infiT 2} 


1 0 


2 0 


(6.9) 
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The exponent v in equation (6.9) is 


v = 2-infi. 


In exactly the same manner, the second integral in equation 
(6.7) can be shown to be 


-(1+infi) infig T 

2 n 2 e X 2n ' 

1+n fi 


where 


2 IT 

T .2 + fi 


'2n 


/ “ | ’■ iV-1 infix, 

1 4>+X 0 1 e dx 


1 + |<|, 1 +X 0 | v - 1 e lnaT l-U 2 +X 0 | v - 1 e inaT 2 . (6. 

Finally, combining the results (6.9) and (6.10) and defining 

1+infi 


10 ) 


*1 = J e 

n ' n 


-infiBn 


2X n (1+infi) 

^ t + r t 

i _ T no 


, . 2 2 In 

1+n fi 


, . 2 n 2 2n ' 

1+n fi 


Equation (6.7) can be expressed as 



fil l r co s nfi (g-0 ) -i 
tt I J rr 1 sin nfi(g-0 n )‘ r 


( 6 . 11 ) 


It is observed that the complex J n is independent of x; for 
each n it is a complex constant dependent on the paramters N, fi , 
and Xq and on the outer limit values of xi and X 2 which are known 
from solving equations (6.5) and (6.6) . In computing the results 
to be discussed later, the integrals in I^ n and l 2 n were evaluated 
numerically using Simpson's Pule for n>_l. I^q and I 20 are easily 
evaluated analytically: 
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t 2 + 27r/n 

I 10 = / (<t>+X 0 ) 2 dT + | (4»J - 4*2 + 2x 0 ^l " <f> 2 )} 

T 1 

= J {N 2 + 2( (I+XqJ^-Xq) 2 } + 2X Q { (I+X^-Xq} (t 1 -t 2 ) 

+ 2XqN { sin ^x^-sin f2x + ^ (cos ^x^-cos ftx 2 }. - D (N) , 


T2+2tt/ 

I 20 = / l+ +x 0 l dT + 

T 1 

= -N { (sin fix^ - sin ftx 2 ) + ^ (cos ^x^ - cos ftx 2 ) } 

+ {(l+x2)>“-x 0 }{if - (T r t 2 )) , 


where the shock relation (6.6) has been used in evaluating I^g . 
It is noted here that, in absence of shocks, i.e., when the ine- 
quality (5.13) is not satisfied, all of the above results which 
determine C n and D n are valid if xj. is simply set equal to x 2 - 


The asymptotic matching of u is now completed by equating 
expressions (6.1) and (6.2) as required by the matching principle 
(ref. 13). This yields the following results which determine com- 
pletely the acoustic field transmitted from the throat region 
into x<0: 


6 = 


A„ = 


A = 
n 


ip = 
n 


e , 



1 

( 2 ) 
' y+l' 

3/2 

2? 


CO 

( 2 ) 
W+i ; 

3/2 

2ira 


0 

n 

1 

&n 

a 2 

2a^ 



oo 


27Ta 


1 n 

2 (Y+l) a 
2 

£ 


j„ -i) . 


n>l , 


(6.12) 
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The results_ (6 . 12) allow complete determination of the com- 
plex constants A^ n in the linear acoustic field (3.11) upstream of 
the throat. The constants Bj_ are then known from equation (4.12) . 
Several important nonlinear effects on the sound field are appar- 
ent immediately. First, the perturbation level in the outer region 
is an order of magnitude smaller than near the throat. This re- 
flects the physically expected transonic trapping of the upstream 
propagating sound. Also, the single harmonic source at the throat 
generates an infinite number of harmonics in the outer region 
which arise, of course, because of progressive nonlinear distortion 
of the waveform as it moves upstream through the near-sonic flow. 
Finally, along with the superharmonics there is generated a steady 
component, or acoustic streaming term, in the outer region which 
is represented by the coefficient Aq in the results (6.12). 

All of the above effects have been shown previously to be 
present even when shock waves are not formed (ref. 5) . However, 
the situation in which the field is continuous is likely to be 
exceptional in practice. The above results demonstrate that in 
general the acoustic field will develop shocks as the throat Mach 
number approaches unity. By far the major physical effect of the 
nonlinearity will be the dissipation of energy between the source 
and the exit associated with these shock waves. In order to eval- 
uate this aspect in detail a discussion of some energy concepts 
will be carried out in the next chapter. 
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7. ENERGY CONSIDERATIONS 


The existence of a shock wave in each period of the inner 
solution gives rise to a dissipation of fluid energy arising from 
viscous and heat-conduction effects modelled by the shock condi- 
tions (3.6). It is of interest to have a convenient measure of 
this dissipation. In linear acoustics, energy is commonly dis- 
cussed in terms of the so-called acoustic energy flux, a quantity 
which is known to satisfy a conservation law involving only first- 
order acoustic quantities. The derivation of this conservation 
law depends upon the fact that the linear acoustic theory is a 
valid leading order approximation in a regular perturbation scheme 
about the basic state. In effect, the common acoustic energy flux 
vector must be the leading order approximation to some physical 
quantity in the exact theory. It is not immediately obvious just 
what this quantity is, and, more importantly, just what is its 
leading order approximation in a singular perturbation scheme such 
as that carried out in the present paper. In the following section 
a consistent representation of perturbation energy is suggested, 
although much of the lengthy detail of its development is omitted. 

A more complete account will appear in a later publication. 


Representation of Perturbation Energy 


The general energy equation for quasi-one dimensional flow 
can be written in terms of the present dimensionless independent 
variables as 


c s (A(x) E (x,t) > t + {A(x)W(x,t) } x = 0, (7.1) 


where E and W are the total fluid energy per unit volume (internal 
plus kinetic) and the total fluid energy flux, respectively; for 
an ideal gas 


E 


P 

Y-l 


+ 


1 

2 


—2 

pu , 


W 


2 




(7.2) 


Now, if the forms (3.3) are substituted into equation (7.1) and 
the basic steady flow equations are used as required, the energy 
equation becomes 

(c ARc?E* } . + {ARUc^W* } = 0 

S 0 t St 


in which 
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E* = y ( y - 1 ' ) ' + IT { ( 1 + P) (1+u) 2 -l} , 

2 

W* = (p+u + pu) + { (l+p) (1+U) 3 -l} . 

Then since ARU, the mass flow rate in the basic flow, and c s are 
constants, the exact energy equation in terms of the perturbations 
p, u, p can be written 



E* + W* 


0 . 


Since, in this work, time periodic solutions p, u, and p are of 
interest it is convenient to time-average this equation over a 
period to obtain 


<W*> x = 0 (7.3) 

in which the notation <> will denote an average over a period. 

The total fluid energy flux, however, is not a very useful 
measure of sound energy, as is well known. In particular, if the 
perturbation quantities are expanded according to equations (3.8) 
it is found that W* contains terms of 0(6), and also that the co- 
efficient of 6 2 contains terms involving P 2 , U 2 , and P 2 which are 
not known if only the leading-order acoustic problem is solved. 

As was sketched in reference 1, and will be discussed in detail in 
a later publication, it proves to be useful to introduce a modified 
energy equation which is of the form 



+ W 

x 


0, 


(7.4) 


in which E=E*-E, and W=W*-W, where E and W can be shown to involve 
contributions to the total energy quantities which are not strict- 
ly attributable to the acoustic perturbation. When the modifica- 
tion is made, equation (7.3) is replaced by 


I = <W> = 0 

x x 


(7.5) 


where 
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w = 


w * - 2 g < 2u+u2 > - Tr=ir (p+u+pu) - MIC ( ITp - 1 ’ 


1 i pu + P (p-p) } + ( (t-Dm -i } 

— t 'k r: lt;u 1+p J T l 5 


(Y-l)G 


(Y-l)G 


pu 


+ {2u 2 + p (3 u 2 +u 3 ) } 


(7.6) 


The quantity W is denoted the dimensionless perturbation 
energy flux and is adopted as the measure of sound energy propaga- 
tion. Its time-average, I, is the dimensionless average power flow 
per period across any station x; it is referred to as the pertur- 
bation intensity. It is emphasized that the conservation law (7.4) 
is an exact relation, which follows directly from the system (3.4) 
as can be verified by constructing suitable combinations of those 
three equations. No expansions of the quantities p, u,and p have 
been used in its derivation. If, however, the outer expansion 
(3.8) is substituted into equation (7.6) it is found that 

2 ° 4 

W = 6 W 2 + 0(6) 

where 


O 



ifpj+U^ (P 1 +M 2 U 1 ) 


G q ( P 0 +U q) ( p o +M u o^ 


(7.7) 


This last expression is exactly the usual acoustic energy flux as 
defined, for example, by Cantrell and Hart (ref. 14). Thus, the 
significance of W is clear: it is precisely that portion of the 

total energy flux in the exact theory whose leading order repre- 
sentation in the outer region is the common acoustic energy flux. 

Physically, W can be shown to be the total fluid energy flux 
less a sum of terms not strictly associated with the propagation 
of sound. These include: (i) the energy flux in the basic steady 

flow, (ii) the flux of internal and kinetic energies resulting 
from convection of the perturbation quantities p, u,and p by the 
steady flow, and (iii) the rate of work done by the perturbation 
pressure against the basic flow velocity and by the basic steady 
pressure against the perturbation velocity. Excluding these con- 
tributions to W has the effect of removing, in equation (7.5) , all 
of the terms in equation (7.3) which are linear in the perturbation 
quantities . 

According to equation (7.5), the perturbation intensity I is 
a constant so long as shocks do not form in the field. When shocks 
appear then | I | must decrease along the propagation direction to 
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account for the energy dissipated by the shocks. In the outer 
region the shocks become simple acoustic discontinuities; it can 
be o shown easily that I, which in the outer region is given by 
<W2>, remains constant in the linearized theory even if these dis- 
continuities are present. 

Thus, in the event that shocks are formed in the inner region, 
the entire effect of these shocks on the sound energy is contained 
in the behavior of W in that region. The inner representation of 
W is complicated to derive because of the dependence on x (and 
thus on X and e) of all the coefficients involving M and G in 
equation (7.6). In view of the fact that W is conserved in the 
shock-free case, it is expected that its leading order inner rep- 
resentation is at order because the results (6.11) have already 
determined that <S=e 2 . This is borne out by introducing the expan- 
sions (4.2) and (4.3) into equation (7.6) and collecting terms in 
the form 


W = I e rn W 1 (X , T) . 

i m 
m=l 

After considerable algebra it is found that, for the present solu- 
tion in which ri=-y 1 and ci=yr 1 , W 1 =0 for m=l ,2,3, while 


W 4 Y+l ^ r 2 +y 2 _y l^ ^ r 2 +y 2 2m l 1J l +Y 2 y l^ * 


It is easily shown that the solution of equation (4.7) for 
r 2 +y 2 is 

r 2 +y 2 = y l + m l y l + 


in which F(T) is an arbitrary function of T whose specific value 
can be determined by asymptotic matching with the outer solution 
but is actually unnecessary for the present purposes. Using this 
result and the inner solution for y^ puts in the form 

M 4 (X ' T > = fyfrr f2<T) - \ (7.8) 

4 

The fact that energy flux at order e in the inner region 
depends on the second order inner solution r2 and U2 is expected. 
It is known (ref. 3) that e^F(T) is equal to the leading order 
inner approximation to the downstream propagating Riemann variable 
for the system (3.4), which is unaffected by the nonlinearity. 
Hence, downstream travelling waves at order 6=e^ i n the outer re- 
gion, arising from continuous reflections from the varying area, 
carry energy at order 6 2 =e^ back into the inner region. 
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Sound Energy Dissipation in the Inner Region 


As discussed above, the change in the sound intensity I 
through the inner region provides a measure of the energy dissi- 
pated by shock waves as the sound field propagates upstream. The 
intensity imparted to the field at the source is calculated to 
leading order using equation (7.8) as 


I 

s 


<WJ(X 0 ,T)> 


2e 


Y+l 


{<f 2 > 


1 

(Y+l) 2 


< ( c(> 2 — Xq- 1 ) 2 > 0 ) » 


where the subscript 0 implies that the average is evaluated at 
X=Xq. On the other hand, the sound intensity in the outer region 
is represented using equation (7.7) by 


out 


£ 

G, 


<( p o+ u o ) (Po +M o u o }> ' 


which can be calculated as an infinite series of harmonics, in- 
cluding a steady streaming term, by use of the real parts of ex- 
pressions (3.11). The difference between I s and I Q ut + s the sound 
power dissipated in each cycle by the shocks occurring between the 
source and the outer region . 


Since I ou t is a constant, and the inner and outer fields have 
already been matched, it follows that I ou -^ must be equal to the 
one-term outer expansion of the one-term inner expression for I, 
which is 


4 i 2e 

lim e <Wt(X,T)> = 

X+-00 * Y 


{<F > - 


(Y+l) 


2 < (<()2_X 0 -1) 2> » } ' 


where the subscript °° indicates that the time average is to be 
evaluated in the outer limit, i.e., using equation (6.3) to relate 
r to T. This determines <F^> in terms of and it follows that 


AI=I -I 


out 


= 2e 4 {< (<j> 2 -X 0 2 -l) 2 > m - < (<J) 2 -X n 2 -l) 2 > n }. (7.9) 


0 


(Y+l) 


The required time averages in the outer limit can be calcula- 
ted as in Chapter 6. For example. 


, 2tt 
T 2 + fl 


, 2ir 
t 2 + 7T 


2 tt 2 v 2 _ C ,}\ 2 _ C \2 2 .. 

ir 0 " x o - - J ^ VdT^ " J ((|) _x o )dT 
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+ {^4 * 0 th-* 2 > } - 
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Use of the shock condition (6.5) puts this in the form 

27 r/ft 

if <^ 2 " X 0 > oo = / (4> 2 -X^)dx -D(N) . 

0 

In exactly the same manner it may be shown that 

2tt/Q t, 

f (<j> 2 -X 2 ) 2 dx - j (cf> 2 -X 2 ) 2 dx 


T 1 


(7.10) 

On the other hand, the time averages at X=Xo can be evaluated 
with reference to figure 3 and to equation (5.25). In the range 
N*<N<N C , which involves two shock branches as depicted in figure 3 
c and d, the integration over the period 0<T<2 ir/ft is broken into 
the intervals 0<T<x+, x + <T<x®, and t^<T<2tt/^. Over the first and 
third of these T=x , while over the middle range 


+ { 


(<j>+ X Q ) 


4X 0 (c(>+X 0 ) 


+ 2X 2 (4>+X q ) 2 } 


2 it . 2 v 2,2 

IT -X 0 > c 


dT , 

2X n <j) ' (x) 
+ „ u 


dx 

<P 2 ( T)“X 2 


and x varies over t_>x>x°. 

Hence 


T + 

T 2 
A** 

X 9 + 2ni/ft 

2tt ^ ,2 2 f ..2 

IT < * “Vo = J ((p " 

X 2 )dx + J (cj) 2 -X 2 ) dx 

+f (4’ 2 -xhdT 

0 

T_ 



+ 2X 0 (<P(t°)-<P(t-) ) . 

The integrals here can be combined into the form 


2 v/ ft 


/ 


( 4> 2- Xq) dx 




(<j> 2 -X 2 )dx + 


/ 


(<j> 2 -X 2 )dx 


L 2 


Then, use of the shock condition (5.29) with D=0 for the lower 
branch of the shock and of equation (6.5) for upper branch yields 
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2tt/Q 


TT < * 2 - x 0 > 0 / “» 2 - x o )dT - D < H > = if < 1> 2 - x § > » • 


It follows in the same manner that 


2 it /S i l l T 

< (<)> 2 -Xq) 2 > 0 - f ($ 2 -X 2 )dT -f (* 2 -x 2 ) 2 dx+/ (* 2 -x 2 ) 2 dr 


- 2X {^ - x 2 *} I 2 

u T 


(7.11) 


Use of the above results in equation (7.9) yields for A I the ex- 
pression 


AI = g4 ^ 3 { < (4> 2 -Xq) 2 > 00 - <(4> 2 -X^) 2 > 0 } (7.12) 

(y+1) it 

in which the averages indicated are given by equations (7.10) and 

(7.11) . It is a simple matter to show that for N<N* equation 

(7.11) retains only the first term on the right side, and that for 
N>N C the terms involving x + and t_ are omitted. 

o o 

In every case, xi, x 2 , X]_, x 2 , x + and x_ are known, and thus 
AI can be simply calculated by substituting equation (5.10) for 
<J>(x). in the numerical results to be discussed later the result 

(7.12) was used in the form 


I out 1 


I , 
out 

to compute the fraction of power input at the source transmitted 
out to the duct exit. 

When no shocks occur then only the first terms on the right 
sides of equations (7.10) and (7.11) appear. In this case, 
equation (7.12) yields AI=0, which reflects the fact that sound 
power is conserved in the absence of shocks. 
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8. DISCUSSION OF RESULTS 


The analytical results derived in the previous chapters have 
been evaluated numerically to study the behavior of sound trans- 
mitted upstream from an acoustic source located at the duct throat, 
Xg=0. All of the following discussion refers to a single duct 
geometry: the area parameter, a, was chosen to be unity, and the 

anechoic exit of the duct was located at x e =-0.5. Of interest are 
the effects on the sound field of variations in the strength of the 
source, the fundamental frequency, and the throat Mach number- In 
the calculations, the acoustic power at the exit, lout/ was com- 
puted by the usual process of time averaging products of real parts 
of the complex solution given in equations (3.11) and (6.12). The 
input power at the source then follows from equation (7.12). 

The total acoustic power I ou t consists of a sum of terms rep- 
resenting the power carried in each harmonic of the outer solu- 
tion. In each case these include a steady streaming component 
n=0. This term actually represents a modification of the basic 
steady flow at the level of magnitude of the perturbation field. 
Hence, a distinction is made in the results between the total 
power transmitted to the duct exit and the fluctuating power, 
which excludes the streaming component and is what would actually 
be measured as sound. It should be noted here that a significant 
streaming effect arises in the outer sound field because of the 
special choice of the source function (eq. (5.10)). Previous work 
in the shock-free case (ref. 4) and preliminary work in the pres- 
ence of shocks indicates that when a simple harmonic source is 
located downstream of the throat region the streaming component 
developed in the field upstream is very small. Of course, in this 
circumstance the perturbation velocity in the throat region is 
markedly different from the form appearing in equation (5.10) . 

One other remark concerning the results is made. The dimen- 
sionless acoustic power actually calculated in the following is 
the actual acoustic power ARUc|<W> divided by a reference power 
corresponding to stagnation conditions, A(0)R s c|. 

Figure 4 is included to show the typical behavior of the shock 
wave in the inner solution. The actual space-time path of the 
shock, T=Hq(X) , in one period is plotted over -20<X£0. The jump 
in ell at the shock, which is a measure of the shock strength, is 
also shown over the same range of X. The shock behaves in typical 
fashion (see ref. 10): its strength grows from zero at X=0 (the 

source location to a maximum at X=«-2, and thereafter decays mono- 
tonically as the wave travels towards the outer region. 

Figures 5 and 6 illustrate the waveform of the acoustic field 
in the outer region for the indicated values of the parameters. 

The wave shapes here are determined by summing 20 harmonics in 
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the series (3.11) for density, pg (fig. 5), and velocity, ug 
(fig. 6) . On each of the figures are shown one period of the wave 
at the exit location x e =-0. 5 , and at the midpoint of the duct, 
x=-0. 25.. The shock wave, calculated using the results of Chapters 
5 and 6 is evident as is the fact that the steady streaming effect 
has grown to non-negligible magnitude. 

On figure 7 is shown the variation in power at the duct exit 
as a fraction of the total sound power introduced at the source 
over a range of values of the source strength parameter N. In this 
case shocks develop in the inner field for N>0. 449. For smaller N, 
total power is conserved. Figure 8 displays the spreading of the 
total power at the exit into harmonic components. The effect of 
energy dissipation in the shocked inner solution is striking. From 
figure 7 at N=6, for example, only about 45% of the source power 
reaches the duct exit in the form of sound. Reference to figure 
8 shows that the power carried in the fundamental harmonic at the 

duct exit is down to less than 30% of that introduced at the source. 

Figure 9 shows the variation of the dimensionless power (re- 
ferred to the reference value discussed previously) with frequency. 
In general it is found that dissipation increases rapidly with 
frequency as is seen here even for the relatively low value of 

N=1.5. This is expected, of course, since the higher the frequen- 

cy the more shock waves occur in a given length of the duct. Fig- 
ure 10 depicts the total and fluctuating power at the exit as frac- 
tions of the input power. At oj= 5 , for example, the combination of 
dissipation and streaming yields a fluctuating power at the duct 
exit only about 20% of that at the source. It is noted, but not 
shown here, that only about 62% of that remains in the fundamental 
harmonic. It should be emphasized that the large attenuations 
seen in the results of figures 7-10 occur even though the throat 
Mach number is only 0.9. Similar results occur for even smaller 
Mach numbers. 

Finally, figures 11 and 12 show the dependence of the outer 
solution on throat Mach number for fixed frequency and fixed source 
strength. It is noted that the effective strength of the velocity 
source at X=0 is eN; its value was fixed here at 0.08. Shocks occur 
in this case for e<0.178 or for M(0) >0.822. On figure 11 are dis- 
played total and fluctuating power at the duct exit as fractions 
of the input power at the source for a frequency w= 2 . The marked 
increase in energy dissipation as the throat Mach number increases 
toward unity is evident. The curves on figure 12 illustrate the 
magnitude of the first three harmonic components of the acoustic 
particle velocity at the duct exit as fractions of the input velo- 
city amplitude eU. As in the results shown in figure 9 these re- 
sults verify the fact that the perturbation quantities at the exit 
actually approach zero as the nonlinearity becomes stronger. 


46 



9. 


CONCLUDING REMARKS 


A nonlinear quasi-one dimensional theory of sound propagation 
through a high-subsonic throat has been developed which includes 
the development of shock waves in the acoustic field. The appli- 
cation of the theory in the present report has been to the special 
case in which the source of sound is located in the region where 
the perturbation field is essentially nonlinear. However, this 
case is of value in understanding the physical propagation process 
in that it illustrates the essential physical phenomena which oc- 
cur as sound interacts with a near-sonic flow, namely the mecha- 
nism by which shocks arise and dissipate sound energy as they prop- 
agate. 

Several conclusions which follow from the preceding analysis 
are worthy of emphasis. These are: 

1. The development of shock waves is the dominant nonlinear 
effect associated with the interaction between sound and flow at 
near-sonic speeds. The shocks give rise to two kinds of physical 
effect on sound. First, they cause a dissipation of acoustic ener- 
gy during the propagation process. Second, they can interact with 
a sound source located in the near-sonic region to cause an acous- 
tic choking of the flow. In this circumstance the source can be 
prevented, during parts of its period, from imparting any sound 
energy to the flow. 

2. In the present application both the nonlinear inner theory 
and the linearized outer theory can be solved analytically. It is 
found from the solution that the shock waves can undergo rather 
complicated motions as they traverse the near-sonic flow region. 

3. The dissipation of sound energy associated with the shocks 
can be discussed in terms of a generalized perturbation energy 
flux expression developed in the report. The usual concept of 
acoustic energy flux is not applicable to the sound field in the 
region of near-sonic flow. It is found that a substantial frac- 
tion of the sound energy introduced by a source can be dissipated 
during the propagation process. 

Current research on the nonlinear behavior of the field from 
a source located well outside the high-subsonic region indicates 
that the nonlinear effects on sound in that case are very similar 
to those in the present study. For example, shock waves arise in 
the near- sonic region in much the same way as in the present anal- 
ysis, and they give rise to energy dissipation of similar magni- 
tude. In addition, the present analytical results prove to be val- 
uable in understanding the numerical analysis required in that 
case. It is hoped that they will also prove to be an aid in in- 
terpreting experimental measurements. 
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APPENDIX 


The exact analytical solution for the Crocco-Tsien duct has 
been discussed in detail by the authors in references 1 and 5. 
The functions Gi n in equation (3.11) can be taken directly from 
section (4.3) of reference 1 or from Appendix D of reference 5 . 
The following parameters are defined: 


3 

£ 

A 


/ T+l x h nw 
K 2 1 K 


sgn x. 


2 + 


2i8 

Y+l ' 


{1 + ^1^ 4 - B 2 
(Y+l) 


/ 


j = (£-l-A)/2, k = (£-l+A) /2 , 
in which K 2 =2a/(Y+1) . Then, if z=(l-K|x|) 2 , 

I-* 1-£ 

G ln = F(-k,-j,2-£;l-z) 


G 2 n = F ( 3 fk,£; 1-z) 

l-£ 

G 3n = ^ 2 T+iB~ (Y+l) P (-k, - j , 2-£; 1-z) 

G 4 n = 2 +i g (Y-l+iB)F(j,k,£;l-z) 

+ (Y+ o } (1-z) F ( j+1, k+1 , £+1 , 1-z) 


where F (a^ , a 2 a^ ; 1-z) is the standard hypergeometric function with 
complex arguments. 

The behavior of the functions G^ n , i=l,...,4, as x-*0 also 
is determined in references 1 and 4 . It is as given by equations 
(3.12) in which the constants are determined from 
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^ 2 ^ a} s 9 n x 



C H = -d + ^ sgnx) 
a 
h 

ino) (Y+l) -2 (Y-l) a 2 sgn x 
inw (y+1) -4a^ sgn x 


2a_) h 


‘y+i j 


sgnx 



REFERENCES 


1. Myers, M. K. ; and Callegari, A. J. : On the Singular Behavior 

of Linear Acoustic Theory in Near-Sonic Duct Flows. Journal 
of Sound and Vibration, vol. 51, no. 4, Apr. 1977, pp. 517- 
531. 

2. Callegari, A. J.; and Myers, M. K. : Nonlinear Effects on 

Sound in Nearly Sonic Duct Flows. AIAA Paper No. 77-1296, 
Oct. 1977. 

3. Myers, M. K. ; and Callegari, A. J. : Transmission of Sound 

Through High-Subsonic Flows in Non-Uniform Ducts. AIAA 
Paper No. 78-1155, July 1978. 

4. Callegari, A. J. ; and Myers, M. K. : Sound Transmission in 

Ducts Containing Nearly Choked Flows. AIAA Paper No. 79-0623, 
March 1979. 

5. Callegari, A. J. : Nonlinear Effects on Sound Propagation 

Through High Subsonic Mach Number Flows in Variable Area 
Ducts. NASA CR- 3138 , May 1979. 

6. Myers, M. K. ; and Callegari, A. J. : Acoustic Shocks in Nearly 

Choked Duct Flows. AIAA Paper No. 80-1019, June 1980. 

7. Nayfeh, A. H. ; Shaker, B. S. ; and Kaiser, J. E. : Computation 

of Nonlinear One Dimensional Waves in Near Sonic Flows. 

AIAA Journal, vol. 16, Nov. 1978, pp. 1154-1159. 

8. Nayfeh, A. H. ; Kelley, J. J. ; and Watson, L. T. : Nonlinear 

Propagation in Near Sonic Flows. Journal of Sound and 
Vibration, vol. 75, no. 3, Apr. 1981, pp. 359-370. 

9. Liepmann, H. ; and Roshko, A.: Elements of Gasdynamics. New 

York, John Wiley and Sons, 1957. 

10. Whitham, G. B. : Linear and Nonlinear Waves. New York, John 

Wiley and Sons, 1974. 

11. Courant, R. ; and Hilbert, D. : Methods of Mathematical 

Physics, Vol. II. New York, Interscience, 1962, Chapter II. 

12. Blackstock, D. T. : Connection Between the Fay and Fubini 

Solutions for Plane Sound Waves of Finite Amplitude. 

Journal of the Acoustical Society of America, vol. 39, 

1966, pp. 1019-1026. 


50 



13 


. Van Dyke, M. : Perturbation Methods in Fluid Mechanics. 

Annotated Edition. Stanford, Parabolic Press, 1975. 

14. Cantrell, R. H.; and Hart, R. W. : Interaction Between Sound 

and Flow in Acoustic Cavities: Mass, Momentum, and Energy 

Considerations. Journal of the Acoustical Society of 
America, vol. 36, 1964, pp. 697-706. 


51 












2 2— 2 O- 1 8- 1 6- 1 4r- 1 2— 10—8 -6 -4 -4 

Axial distance, X 


Figure 


4. -Shock strength and path: N-5, £-0.1, w 2 
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Figure 5. -Time history of perturbation 
density Po : N=2, w=2, e=0.1, 

x =0.5, 20 harmonics, 
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Figure 6. -Time history of perturbation 

velocity Uq : N=2, to=2 , e=0.1, 

x =0.5, 20 narmonics. 
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Figure 7. -Fraction of input power reaching exit 
vs. source strength: e=0.1, w=2. 
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Figure 8. -Fraction of total power at exit in harmonics 
vs. source strength: e=0.1, w=2. 
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Figure 11. -Fraction of input power reaching exit 
vs. throat Mach number: eN=0.08, to=2. 
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